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ABSTRACT 

The thermal diffusivity and thermal conductivity for a crushed 
olivine basalt were determined from transient state data. Values were 
obtained over a temperature range of — 100 to 200 °C in vacuums of 
5 X 10 3 and 5 X 10 _G mm Hg as well as at atmospheric pressure. 
A — 150 mesh material at a density of 1.14 g/cm 8 had a thermal con- 
ductivity of 3.9 X 10~ G cal/cm sec °C at 100 °C when measured in a 
vacuum of 5 X 10~ (i mm Hg. This was approximately one hundred 
times lower than the values obtained for the same material measured 
at atmospheric pressure. Increasing the density to 1.57 g/cm 3 increased 
the thermal conductivity by approximately 60% in both air and 
vacuum. Over the range studied, the test temperature had very little 
effect on thermal conductivity in air but showed more of an effect 
when the material was placed in a vacuum. 


I. INTRODUCTION 


In preparation for manned landings on the Moon, re- 
search aimed at determining the physical and mechanical 
properties of the postulated lunar surface is now in 
progress. There have been extensive astronomical studies 
of the Moon and its surface; on the basis of these, 
it is proposed that the Moon is composed of igneous 
rocks similar to certain types that are found on Earth. 
Some measurements (Ref. 1, 2, 3), however, indicate 
that the thermal conductivity of the lunar surface is 
extremely low — much lower than that of any known 


types of consolidated rock. To account for this, it is 
postulated that large portions of the lunar surface are 
covered with a highly porous material, probably powdery, 
perhaps lightly sintered. 

The investigation described here is one part of a pro- 
gram designed to evaluate the properties of a powdered 
rock simulating the postulated lunar surface material. Its 
purpose was to measure the effects of vacuum on the 
thermal diffusivity and conductivity of rock powder. 
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II. MATERIALS TESTED 


The material used in the initial runs, which were de- 
signed to check equipment operation, was a 70-mesh 
commercial grade 99% silica foundry sand. The particle 
size distribution of this material is given in Table 1. The 
moisture content was less than 0.03%. 

The material selected to simulate lunar material for 
this study was an olivine basalt collected from Pisgah 
Crater, San Bernardino County, California. The large 
volcanic bombs were crushed to about 3-in. size using a 
hydraulic press. Any pieces which appeared to be con- 
taminated or which showed signs of weathering were 
discarded. The remaining material was passed through 
a Gates jaw crusher and a stainless-steel hammer mill- 
type pulverizer to reduce it all to —35 mesh material. A 


detailed mineralogical description of this material and 
the crushing procedure is given in Ref. 4. 

Most of the —35 mesh material was screened into more 
closely sized fractions; nominally —35 +48, —48 +65, 
— 60 +100, —100 +150, and —150 mesh. Actual screen 
analyses of the fractions used for the thermal conduc- 
tivity tests are shown in Table 1. 

The “as received” density of the rock as measured by 
water displacement was 2.83 g/cm 3 at 24° C. The density 
of the crushed material was 2.97 g/cm 3 . All of the 
screened fractions were stored in closed containers. 
Measurements made during the program showed that 
actual moisture content varied from day to day but 
was always less than 0.1%. 


Table 1 . Particle size distribution of materials tested 


Tyler 

mesh size 

Particle 
size, ju 

99% silica 
sand 

% 

Olivine basalt, % 

Nominal 
mesh size 
-35 

Nominal 
mesh size 
-35 + 48 

Nominal 
mesh size 
-150 

+ 28 

>590 

0.00 

0.05 

0.02 

0.00 

-28 + 35 

<590>420 

0.11 

0.27 

0.59 

0.02 

-35 + 48 

<420>297 

9.11 

8.05 

67.65 

0.03 

-48 + 65 

<297>21 0 

21.12 

10.78 

28.80 

0.03 

-65 + 100 

<21 0> 1 49 

37.45 

10.26 

1.01 

2.15 

-100 + 150 

< 1 49> 1 05 

25.75 

9.66 

0,21 

2.85 

-150 + 200 

<105 >74 

5.61 

9.68 

0.10 

7.98 

-200 + 270 

<74 >53 

0.37 

11.21 

0.09 

19.47 

-270 + 325 

<53 >44 

0.17 

4.94 

0.04 

14.03 


<44 >35 


3.98 

0.05 

4.13 


<35 >25 


9.11 

0.04 

10.70 


<25 >12 


10.88 

0.26 

17.25 


<12 >7 


6.31 

0.29 

9.95 


<7 >3 


3.88 

0.20 

8.82 


<3 


0.61 

0.03 

2.08 

Total 


99.69 

99.87 

99.28 

99.49 

Loss in screening 


0.31 

0.13 

0.72 

0.51 



100.00 

1 00.00 

100.00 

1 00.00 
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III. TEST APPARATUS AND EXPERIMENTAL PROCEDURE 


It was expected that very low thermal conductivity 
values would be obtained for the loosely packed mate- 
rials in vacuum; thus, conventional equipment could not 
be used for these measurements, and a special apparatus 
was designed and built. 

A general view of the equipment used is presented in 
Fig. 1; details of the specimen holder configuration are 
shown in Fig. 2. This apparatus was designed to measure 
the rate of change of the temperature at points along a 
diameter of a cylindrical specimen caused by a sudden 
stepped change in the outer surface temperature. To 
minimize errors due to heat conduction along the tem- 
perature sensors, very thin (0.003-in. D) unshielded 
chromel-alumel thermocouples were strung lengthwise 
through the specimen holder (Fig. 2). One couple was 
located on the centerline, two at midradius, and two at 



HEATING 

COIL 


COOLING 

COIL 


SPECIMEN 

THERMOCOUPLES 


Fig. 1. Thermal diffusivity test equipment 


the outer surface of the specimen, which was also the 
inner wall of the cylindrical specimen holder. A 1500-w 
hermetically sealed heating coil and a copper cooling coil 
were interwound on the outside of the specimen holder. 
A two-piece polytetrafluoroethylene closure was pro- 
vided for the bottom of the holder. To increase thermal 
resistance between the thermocouples and the copper 
specimen holder, the bottom closure was so designed 
that it made only a line contact with the copper holder 
and also with the thermocouples. 

A similar closure was used as a top cover during some 
of the preliminary runs with the thought that it would 
reduce radiation losses from the upper surface of the 
test specimen. However, it was determined that the meas- 
ured diffusivities were not affected by the presence of 
such a cover; so, during all subsequent runs, the cover 
was not used, thus allowing visual observation of the 
surface of the specimen during the entire test. 

The holder assembly and the supporting structure were 
attached to a double-flanged spacer which contained the 
vacuum leadthroughs for the thermocouples, heater 
cables and liquid-nitrogen supply lines. The vacuum 
system consisted of a 12-in. glass jar, a liquid-nitrogen- 
trapped 4 in. oil diffusion pump (210 1/sec pumping 
speed at 8 X 10~ 4 mm Hg), and a 13.3-cfm mechanical 
pump. Pressures on the order of 5 X 10 -6 mm Hg could 
be reached in about 2 hr and maintained for the dura- 
tion of the test run. 

Earlier studies of a similar material had indicated that 
the packing density had a very significant effect upon the 
properties ( Ref. 4 ) . Thus, it was necessary to use extreme 
care in loading test material into the specimen holder so 
that packing density could be known and controlled. To 
accomplish this, a known weight of the test material was 
loaded into the specimen holder for each run. When it 
was desired to conduct the experiment with a low pack- 
ing density, the material was carefully spooned into the 
holder and every precaution taken to prevent the holder 
from being bumped or vibrated. The height of the mate- 
rial in the holder was recorded and observed during the 
entire test to assure that vacuum-pump vibration or other 
external causes did not produce a significant change in 
the packing density. When it was desired to conduct the 
experiment with a high packing density, the material was 
loaded in the holder, which was then tapped with a soft 
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hammer until no further decrease in height of the mate- 
rial in the holder could be noted. Again this height was 
recorded and observed during the entire test to assure 
that degassing during evacuation did not significantly 
change the surface level of the test material. 

After loading, the temperature of each thermocouple 
was read to ensure that the entire unit was thermally 
stabilized. When the test was to be run in vacuum, the 
assembly was first evacuated to the desired pressure 
range. The temperature of the specimen holder was then 
changed suddenly by energizing the heater or flowing 
liquid nitrogen through the cooling coil as required for 
a change of approximately 75° F. This sudden tempera- 
ture change was performed manually to assure maximum 
rate of change. An automatic controller was then used to 
maintain the outside wall of the specimen at the new 
temperature level. The outputs of the thermocouples in 


the test specimen (Fig. 2) were recorded on a multipoint 
recorder. 

Each test was continued until the temperature at the 
center of the specimen was within about 10°F of 
the wall temperature, at which time the temperature of 
the holder could again be changed suddenly to extend the 
run to a new temperature range. In air, less than 1 hr 
was required for the center temperature to reach the 
desired point, whereas in vacuum the time required was 
on the order of 36 hr. Specimen temperature vs time 
plots were prepared for each test. A typical example of 
a plot for a test in vacuum is shown in Fig. 3. In this test 
the initial and final temperature excesses (T 2 — T r and 
X 2 — T; Fig. 4 and Table 2) measured for the center 
thermocouple were 92 and 9°F, respectively. For the 
mid-radius thermocouple, the values measured were 91 
and 5°F, respectively. 
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Rg. 3. Temperature vs time for vacuum test 9-1 1 
(olivine basalt —35 +48 mesh at a 
density of 97 Ib/ft 3 ) 




0.80 


0.60 


0.40 


0 . 20 | 


0.10 


0.08 


0.06 



FOR TEMPERATURES MEASURED AT: 
CENTER (REF. 5) 

MJDRADIUS (REF. 6) 


U= RATIO OF TEMPERATURE 
EXCESSES, i.e. 1 (r 2 ^)/(72-r l } l 
WHERE r- CENTER OR MIDRADIUS 
TEMPERATURE AT TIME f t 
TEMPERATURE OF SPECIMEN 
BEFORE THERMAL PULSE, AND 
7^- SPECIMEN SURFACE 
TEMPERATURE AFTER THERMAL 
PULSE; 5 = SPECIMEN RADIUS 
AND a = THERMAL DIFFUSIVITY 


0.8 


af/s 1 


Fig. 4. Change in temperature of infinite cylinder 
caused by sudden change of surface 
temperature 
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IV. EVALUATION OF MEASUREMENTS 


Methods for calculating the thermal diffusivity of a 
material from such transient state data are described in 
several texts. Jakob (Ref. 5) presents a graphical method 
for determining the diffusivity when the (T 2 — T)/(T 2 —T X ) 
ratio ( 17) of the center of a specimen of given geometry 
is known at any particular time. Russell (Ref. 6) de- 
scribes this solution for the midradius position. These 
graphical methods for an infinitely long cylindrical sample 
are illustrated in Fig. 4. It was found from the error 
analysis ( see Appendix ) that the most accurate values of 
diffusivity are obtained if calculated from the segment 
of the time-temperature curves (Fig. 3) where U is in 
the range 0.2 to 0.5. Four diffusivity values were calcu- 
lated for each test run; that is, two values from the center 
temperature curve at the points where U was near 0.2 
and 0.5 and two from the midradius temperature curve 
at the same U values. The four diffusivities so obtained 
were averaged, and these average values are listed in 
Tables 2 through 5, inclusive. 

Thermal conductivity can be calculated from diffusivity 
values if the density and specific heat are known. The 
packing density was measured for each test as indicated 
previously. Consistent specific-heat values for pure silica 
over a wide temperature range are available in the litera- 
ture (Ref. 7, 8). These data are shown in Fig. 5, and 
were used to make the calibration calculations. 



Fig. 5. Specific heat vs temperature (silica sand 
data from Ref. 7 and 8) 


Specific-heat data for olivine basalt (or similar mate- 
rials) are limited, extend over a narrow temperature 
range, and may be sensitive to the composition. Because 


Table 3. Thermal diffusivity and conductivity of —35 +48 mesh olivine basalt 0 


Test 

No. 

Tost 

pressure 
mm Hg 

Test temperatures, °F 

Density 

Specific 

heat 

Btu/lb °F 
cal/g °C 

Thermal diffusivity 

Thermal conductivity 

T i 




g/cm 3 

ft 2 /hr 

cm 2 /sec 

Bfu/hr ft °F 

cal/sec cm °C 

BEI 

760 


160 

115 

85 


0.230 

5.38 X IO - 3 

1.39 X 10-3 

1.05 X 10~* 

4.34 X 10~ 4 

ESS 

760 


239 

200 

85 


0.240 

5.12 X 10~ J 

1.32 X 10~3 

1.04 X 10-' 

4.30 X lO" 4 

9-1 

760 


160 

115 

97 


0.230 

5.21 X TO- 3 

1.34 X 10-3 

1.16 X lO' 1 

4.80 X 10~ 4 

9-1 A 

760 


240 

200 

97 


0.240 

4.83 X 10-3 

1.24 X TO" 3 

1.12 X 10-* 

4.63 X 10- 4 

10-10 

760 


-32 

18 

85 


0.215 

4.96 X 10-3 

1.28 X 10-3 

9.06 X lO' 2 

3.74 X 10“ 4 

1 0-1 0A 

760 

—31 

-132 

-82 

85 


0.193 

4.84 X 10-3 

1.25 X 10~3 

7.94 X lO' 2 

3.28 X 10~ 4 

10-9 

760 

66 

-30 

18 

97 


0.215 

5.26 X IQ' 3 

1.35 X 10-3 

1.09 X 10-i 

4.50 X 1 0~ 4 

10-9A 

760 

-31 

-132 

-82 

97 

1.56 

0.193 

5.07 X 10-3 

1.31 X 10-3 

9.49 X lO" 2 

3.92 X 10~ 4 

9-11 

5 X 10^ 

70 

162 

116 

97 

1.56 

0.230 

1.32 X 10~ 4 

3.40 X 1 0~3 

2.94 X lO' 3 

1.21 X lO' 5 

9-12 

5 X 10-3 

154 

240 

197 

97 

1.56 

0.240 

1.76 X 10~ 4 

4.54 X lO- 5 

4.10 X lO' 3 

1.69 X lO- 5 


5 X 10-* 

69 

161 

115 

97 

1.56 

0.230 

1.15 X 10~ 4 

2.96 X 10-3 

2.56 X 10“3 

1.06 X lO' 3 


5 X 10-* 

148 

238 

193 

97 

1.56 

0.240 

1.36 X 10~ 4 

3.50 X 10-3 

3.16 X lO" 3 

1.30 X lO" 5 


5 X 10-* 

73 

-32 

20 

97 

1.56 

0.215 

1.11 X 10" 4 

2.86 X 10“3 

2.31 X lO" 3 

9.55 X 1 0~ * 


5 X 10-* 

-14 

-128 

-71 

97 

1 

1.56 

0.195 

1.52 X I0~ 4 

3.92 X IQ” 5 

2,87 X 10-3 

1.18 X lO” 5 


•See footnote* to Table 2. 
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of this lack of reliable data, determinations of the specific 
heat of the olivine basalt used were made between room 
temperature and 100, 200, and 300° F. These determina- 
tions were made using a calorimetric method in accord- 
ance with American Society for Testing and Materials 
Standard Procedure C351-61. These measured specific 
heat values are plotted in Fig. 5. It has been reported 
( Ref. 9 ) that within limits, the variation of specific heat 
with temperature is similar for various types of crushed 
rock. With this in mind, the curve through the measured 
points for the olivine basalt (Fig. 5) was extrapolated to 
give specific-heat values over the temperature range used 
in this investigation. Thermal conductivity values listed 
in Tables 2 through 5 were calculated from the average 
thermal diffusivity values determined graphically (Fig. 
4) and the measured densities and specific heat values 
(Fig. 5). 

For a few of the tests, thermal conductivity was also 
calculated from the experimental data using a digital 
computer. The computer solution avoids the assumption 
made in the graphical method that conductivity and spe- 
cific heat are constant over the temperature range of a 


single test. Values obtained with the computer were in 
agreement with those obtained graphically. 

In order to facilitate analysis of the test data, the 
diffusivity and conductivity values are shown graphically 
in Fig. 6 through 10, inclusive. The last plot (Fig. 10) 
shows the thermal diffusivity and conductivity for the 
three lots of basalt measured at approximately the same 
packing density. The effect of particle size may be readily 
evaluated. 

An error analysis of the measuring and calculation pro- 
cedure was also made ( see Appendix ) . The chief sources 
of error were found to be heat flow through the ends of 
the test specimen, lack of complete temperature uni- 
formity at the start of a test, and uncertainties in the 
specific heat. Less important sources were temperature, 
time, and position measurements, heat flow along thermo- 
couples and bottom closure, and variation of thermal 
properties with temperature during the course of a test. 
The estimated error from all of these sources in each 
value of diffusivity is about 10%, in each value of con- 
ductivity about 15%. 


Table 4. Thermal diffusivity and conductivity of -35 mesh olivine basalt* 2 


Test 

Test 

pressure 
mm Hg 

Test temperatures, °F 

Density 

Specific 

heat 

Thermal diffusivity 

Thermal cdnductivity 

No. 

r , 

T 2 

Lvg 

Ib/ft 5 

g/cm 3 

Btu/lb °F 
cal/g °C 

ftVhr 

cm 2 /sec 

Btu/hr ft°F 

cal/sec cm°C 

4-25 

760 

73 

124 

99 

93 

1.49 

0.227 






4 - 25 A 

760 

122 

176 

149 

93 

1.49 

0.234 

4.16 X 10~ 3 

1.07 X 

IO- 3 

9.05 X 10~ 2 

3.74 X 10~ 4 

4 - 25 B 

760 

174 

228 

201 

93 

1.49 

0.240 

5.91 X 10- 3 

1.53 X 

10- 3 

1.32 X 10 "' 

5.45 X 10“ 4 

8-28 

760 

89 

161 

125 

113 

1.65 

0.231 

6.12 X 10~ 3 

1.58 X 

10~ 3 

1.59 X 10 -' 

6.57 X 10~ 4 

8 - 28 A 

760 

162 

242 

202 

113 

1,65 

0.240 

6.45 X 10~ 3 

1.66 X 

10~ 3 

1.75 X 10 -' 

7.23 X 10~ 4 

4-27 

760 

65 

23 

44 

93 

1.49 

0.219 

6.97 X 1 0~ 3 

1.80 X 

10~ 3 

1.42 X 10 -' 

5.87 X 10~ 4 

4 - 27 A 

760 

12 

-32 

-10 

93 

1.49 

0.209 

7.55 X 1 0~ 3 

1.95 X 

10~ 3 

1.47 X 10 -' 

6.07 X 10~ 4 

10-13 

760 

73 

-31 

21 

93 

1.49 

0.215 

4.76 X 10~ 3 

1.23 X 

10~ 3 

9.52 X 10-2 

3.94 X 10~ 4 

10 - 13 A 

760 

-34 

-131 

-83 

93 

1.49 

0.193 

4.57 X 1 0~ 3 

1.18 X 

10~ 3 

8.20 X 10-2 

3.39 X 10~ 4 

10-16 

760 

70 

-31 

19 

122 

1.95 

0.214 

5.67 X 1 0 3 

1.46 X 

10~ 3 

1.48 X 10 -' 

6.12 X 10~ 4 

1 0-1 6 A 

760 

-39 

-132 

-86 

122 

1.95 

0.192 

5.50 X 1 0~ 3 

1.42 X 

10- 3 

1.29 X 10 -' 

5.33 X 10~ 4 

7-1 1 

5 X 10- 6 

73 

159 

116 

93 

1.49 

0.230 

6.49 X 1 0~ 5 

1.67 X 

10~ 3 

1.39 X 10“ 3 

5.75 X 10 -* 

7-13 

5 X 1 0 -6 

146 

240 

193 

93 

1.49 

0.239 

7.63 X 10~ 3 

1.97 X 

10~ 3 

1.70 X 10- 3 

7.03 X 10 ~* 

7-17 

5 X 1 0 -6 

71 

167 

119 

113 

1.65 

0.230 

9.93 X 1 0 -5 

2.56 X 

10- 3 

2.58 X 10- 3 

1.07 X 10” 3 

7-19 

5 X 10 ^ 

160 

241 

201 

113 

1.65 

0.240 

7.77 X 1 0 -5 

2.00 X 

10~ 3 

2.10 X 10“ 3 

8.68 X 10" 6 

7-26 

5 X 10~ 4 

71 

156 

114 

113 

1.65 

0.229 

8.92 X 10 - J 

2.30 X 

10“ 3 

2.31 X 10- 3 

9.55 X 10 ~* 

7-27 1 

5 X 1 0 ~* 

140 

236 

188 

113 

1.65 

0.239 

9.51 X 10" 3 

2.45 X 

10~ 3 

2.57 X 10- 3 

1.06 X 10~ 3 

10-24 : 

5 X 1 0 “* 

66 

-30 

18 

93 

1.49 

0.214 

5.00 X 10- 3 

1.29 X 

10~ 3 

9.95 X 10 ~< 

4.11 X 10“ 6 

11-1 

5 X 10" 4 

75 

-124 

-24 

109 

1.75 

0.206 

4.56 X 10~ 3 

1.17 X 

10" 3 

1.02 X 10“ 3 

4.22 X 10 ~* 

12-18 

5 X 1 0 “* 

72 

-32 

20 

98 

1.57 

0.215 

5.00 X 10“ 5 

1.29 X 

10~ 3 

1.05 X 10- 3 

4.34 X 10 ~* 

12-21 

5 X 10- 4 

-15 

-120 

-68 

98 

1.57 

0.196 

4.71 X 10“ 3 

1.22 X 

10- 3 

9.05 X lO " 4 

3.74 X 10 -* 


a See footnotes to Table 2 . 
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V. DISCUSSION OF RESULTS 


The thermal properties of silica sand determined in air 
(Table 2 and Fig. 6) are in good agreement with litera- 
ture values (Ref. 10). At reduced pressures (10~ 3 and 
10 -6 mm Hg range), the measured diffusivities and con- 
ductivities were roughly 100 times lower. There were no 
comparable thermal property data at these vacuum levels 
found in the literature. Tests in poorer vacuum environ- 
ments are reported (Ref. 11 ), and it is interesting to note 
that the thermal conductivity of sand at 10" 1 and 10 -2 
mm Hg is not greatly different from the values obtained 
in the present investigation at 10“ 3 and lO^ 6 mm Hg. The 
major change in thermal properties occurs when die 
interstitial air pressure drops from atmospheric down to 
about 10~ 2 mm Hg; beyond this, the thermal diffusivity 
and conductivity are relatively constant. 

The results of the tests on the olivine basalt are given 
in Tables 3 through 5 and Fig. 7 through 10. Again the 
thermal diffusivity and conductivity as measured in vac- 
uum were approximately 100 times lower than the values 
measured in air. These data also show that the effects of 
vacuum are very much greater than the effects of tem- 
perature, packing density, or particle size of the test 
material. 

It is interesting to note that the values of thermal con- 
ductivity obtained for the olivine basalt in vacuum are in 
good agreement with the value (3 to 10 X 10" 6 cal/cm 
sec °C) of thermal conductivity calculated for the lunar 
surface from measurements made during a lunar eclipse 
(Ref. 1,2). 

The effects of the vacuum may be explained by con- 
sidering the thermal conduction process. For the tests in 
air, the conduction is through the particles of basalt and 
through the air which fills the voids. For the tests in 
vacuum, there is conduction through the basalt particles 
but very little through the evacuated voids. Thus, a major 
portion of the thermal energy must be transferred through 
the voids by radiation, which is very low at the tempera- 
ture used in this study. Therefore, the thermal conduc- 
tivity as measured in vacuum would be lower than the 
value measured in air. 

The degree of vacuum in which the measurement was 
made would affect the amount of energy transferred by 


conduction through the voids, but this would be small 
except for very poor vacuums. While the equipment used 
was not particularly suited to running tests at various 
levels of vacuum, it was possible to run some tests at two 
levels. The results of these (Tables 2 and 3) show that 
increasing the pressure from 10~* to 10~ 3 mm Hg increased 
the conductivity by a small amount. It would therefore 
seem reasonable to assume that at pressures below KT 6 
mm Hg, the thermal conductivity will not be markedly 
affected by the pressure. 

The effect of test temperature on the values obtained 
during the tests in air is small (Fig. 7-10). This is what 
would be expected, as the thermal conductivity of the 
basalt particles and of the air in the voids between the 
particles would not change markedly over the range of 
temperatures studied. Since the test temperature would 
affect the radiation through the evacuated voids, it should 
affect the thermal conductivity as measured in vacuum. 
The lower the test temperature, the lower the conduc- 
tivity would be. Although the test data are limited, they 
do seem to show some effect of temperature (Fig. 7-10). 

For a given particle size distribution, the higher the 
packing density, the lower the percentage of voids and 
the more continuous the conduction path. Thus, it would 
be expected that the thermal conductivity in both air and 
vacuum should increase with an increase in packing den- 
sity. Again, the available data are limited; however, they 
do indicate this behavior. 

The effects of the particle size distribution are not 
as clearly understood. The data available from the tests 
run during this investigation are shown in Fig. 10. From 
these, it would appear that the particle size distribu- 
tion has very little effect upon the thermal conductivity 
as measured in air but does have an effect on the value 
measured in vacuum: material with a sharp and coarse 
distribution ( — 35 -f 48 mesh) has a higher conductivity 
than the material with a broader distribution of particle 
sizes. Because of other factors which could also be 
affecting these data, they should not be considered 
as conclusive. 
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VI. CONCLUSIONS 


From the results obtained during this investigation, the 
following conclusion may be stated: 

1. The thermal conductivity of the crushed olivine 
basalt used was approximately 100 times lower 
when measured in vacuum than when measured in 
air at atmospheric pressure. The value measured 
in vacuum is in good agreement with the value for 
the lunar surface calculated from astronomical data. 

2. Increasing the pressure from 5 X 10^ 6 to 5 X 10 -3 
mm Hg had no marked effect on the thermal con- 
ductivity of the crushed basalt. 

3. For the —150 mesh material the thermal conduc- 
tivity in air and in vacuum was increased approxi- 


mately 60% at all test temperatures when the 
packing density was increased from 1.14 to 1.57 
g/cm 3 . 

4. Decreasing the average temperature of the crushed 
basalt specimen from 100 to — 70°C caused a de- 
crease in the thermal conductivity. 

5. For the particular distributions used, the particle 
size had a greater effect on the values of thermal 
conductivity measured in vacuum than on the 
values measured in air. 

6. The thermal conductivities of the olivine basalt and 
the silica sand are not markedly different. From 
this, it is concluded that composition of the crushed 
material has only a minor effect on its thermal 
conductivity. 


14 



JPL TECHNICAL REPORT NO. 32-368 


APPENDIX 

Error Analysis 


I. SUMMARY 


An analysis was made of the various sources which 
introduce errors into the calculated values of thermal dif- 
fusivity and conductivity. This indicated that the over-all 
accuracy of the individual diffusivity values was probably 
about 10% and of the conductivity values calculated 
manually about 15%. For conductivity values calculated 
by machine, the accuracy was probably about 10% . 

Important in the error analysis is the fractional 
temperature 

v = T^T, (A1 > 


where T = temperature at center or midradius and T x 
and T 2 are as defined in Table 2. To keep the errors in 
the calculated properties low, only values of U between 
0.2 and 0.5 were used for manual calculations and only 
values between 0.2 and 0.9 for machine calculations. 

The amounts of error believed to be contributed by 
various sources are tabulated in Table A-l. These num- 
bers are combined by taking the square root of the sum 
of their squares. 


II. DETAILS 


A. Errors Arising From Time Measurements 

Diffusivities ( a ) are calculated from the quantity at/s 2 
of Fig. 4, where t — time from start of test and s = 
cylinder radius. An error in time measurements. At, then 
produces a corresponding error in diffusivity, Aa, given 
by 



The accuracy of time measurements was about 10 sec. 
The time of run t was typically 16, 24, and 30 min for 
tests in air at U = 0.5, 0.2, and 0.1, respectively; in 
vacuum, corresponding typical times were 12, 18, and 
25 hr. At U = 0.9, in air, a typical center time was 12 
min, a midradius time 6 min; in vacuum, a typical time 
was 4V2 hr. Substituting these values in Eq. (A-2) gave 
100 Aa/a, the percent error in a listed in Table A-l. 


where p = density of sample and c — specific heat. 

= (a-4) 

k a 

or the percent error in conductivity is equal to the per- 
cent error in diffusivity. 

In both manual and machine computations, the values 
obtained for k and a are based on several experimental 
measurements: temperature readings at several times 
and at both center and midradius. Because of the aver- 
aging involved in this, the error in the resulting thermal 
quantity will be less than the error in an individual 
reading by a factor of 1 \fn, where n = number of meas- 
urements. 


As the thermal conductivity k is related to diffusivity 
by 


B. Errors Arising From Temperature Measurements 

The slope of the curves (Fig. 4) indicates that an error 
in fractional temperature. At/, produces a corresponding 
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fractional error in a t/s 2 , which is given by [ — A (at/s 2 )/ 
(a t/s 2 ) ]/A(7 = 4.1, 2.6, 2.1, 2.7, and 4.2 for U = 0.9, 0.8, 
0.5, 0.2, and 0.1, respectively. As 

ia = 

a a t/s 2 

these ratios give the fractional error in diffusivity and, by 
Eq. (A-4), in conductivity, resulting from an error in 
temperature measurement. Temperature differences were 
measured to 1°F; this includes errors in thermocouples 
and recorder calibration. The difference T 2 — T was 
typically 50° F, so A C7, by Eq. (A-l), was about 2%. 
Corresponding errors in a are listed in Table A-l. Again, 
these estimates could be reduced by l/\/n. 


C. Errors Arising From Position Measurements 

Errors in radial dimensions As, such as in the thermo- 
couple positions, contribute to errors in a in two ways. 
First, they mean that the applicable curve (Fig. 4) is not 
that for center or midradius but a curve for some slightly 
different position. The Figure suggests that the corre- 
sponding ( Aa/a ) /( A s/s ) = [A(a t/s 2 )/(a t/s 2 )]/ 
(A s/s) is about — 1, —0.6, and —0.4 at U = 0.5, 0.2, 
and 0.1, respectively. At U — 0.9, it is about —1 at center 
and —2 at midradius. Second, for any given value of 
a t/s 2 , 

±JL=*A± (A-6) 

a s 

Summing these effects algebraically, one finds that for 
midradius readings at 17 = 0.9, the two error contribu- 
tions cancel to a first approximation. For center readings 
at U = 0.9, (Aa/a) /(As/s) is about 1; for readings at U — 
0.5, 0.2, and 0.1, this ratio is, respectively, 1, 1.4, and 1.6. 

The thermocouple positions are believed to be accurate 
to 0.002 in.; the cylinder dimensions are more accurate. 
Dimension As may then be taken as 0.002 in.; the nominal 
radius s = 1 in. The value A s/s = 0.002 was used as the 
basis for the 100 Aa/a values of Table A-l. 


D. Errors Arising From Heat Flow Through Ends 

The equations and curves used in calculating thermal 
diffusivity and conductivity from the experimental data 
assume that the sample cylinder may be considered in- 
finitely long: that axial heat transfer does not affect the 
temperatures at the midlength. The error due to the finite 
length of the cylinder, which permits heat flow through 
the ends, can be estimated by using the superposition 
principle. One may consider the resulting change in 


temperature as that occurring in a superposed infinite flat 
plate placed perpendicular to the cylinder axis, with its 
midplane at the cylinder midlength and its surfaces at 
the cylinder ends. End conditions on the cylinder cor- 
respond to surface conditions on the plate. Tables are 
available (Ref. 6) showing the variations of the fractional 
temperature at the plate midplane, U p ut «, as a function of 
the quantity a t/s 2 pU te, for an abrupt change of tempera- 
ture at the plate surfaces. Here, s P i at e is half the plate 
thickness. For the dimensions used in the experiment. 


s pute = 2 in. =2 s 
At any time, then, 


at _ Va at 

S* Plate ~ S 2 


(A-7) 

(A-8) 


From Fig. 4, one may read the values of a t/s 2 cor- 
responding to various values of U at cylinder center and 
midradius and, by Eq. (A-8), get the corresponding 
values of a t/s 2 P i a te. Reference 6 then gives corresponding 
values of U pl&te . The quantity 1 — U piiLte may be consid- 
ered as the change in U that would result if the cylinder 
ends were at all times at the temperature of the cylindri- 
cal surface. 


In practice, it seems likely that the end temperature 
during a test remained closer to the initial temperature 
than to the cylindrical surface temperature; tests were 
always made by stepping the cylindrical surface tempera- 
ture away from room temperature, and the ends could 
“see’* (exchange radiation with) mostly objects near room 
temperature. The fractional temperature difference be- 
tween end and midlength would then be near zero at the 
start of a test and 1 — 17 when a reading was taken, or 
(1 — U)f 2 on the average. For ends at the cylindrical 
surface temperature, this difference would be 1 at the 
start and 17 when a reading was taken, or (1 + U )/ 2 on 
the average. The value of 1 — 17 pIate , calculated as above, 
was therefore multiplied by (1 — 17)/(1 -h 17) to give 
an estimate of the error in U . 

The values obtained for Al7 = (1 — I7 pl a te ) (1 — 17)/ 
( 1 4* U ) were 0.00002, 0.001, 0.06, and 0.023 at center and 
0.000005, 0.0002, 0.04, and 0.013 at midradius for 17 = 
0.9, 0.5, 0.2, and 0.1, respectively. These values were 
multiplied by the ratios of — (Aa/a) /A 17 stated in 
Section II-B to give fractional errors in diffusivity; the 
errors are listed, as percents, in Table A-l. 


E. Errors Arising From Heat Flow Through 
Thermocouples 

Heat flow along the thermocouple wires will change 
the temperature at the junctions and introduce a tem- 
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perature error. Jakob (Ref. 5) gives a solution for the 
steady-state temperature distribution of a wire along the 
axis of a finite insulating cylinder whose cylindrical sur- 
face is at constant temperature. Inserting the values 
pertinent to this experiment, the solution is 

— T 0 = 2(T 0 — T i) sinh 2 ^ Opiate yjj^r 2 lny 

(A-9) 


where 

T l = temperature of wire at end of cylinder 
T 0 = temperature of wire at midlength of cylinder 
Opiate = half the axial length = 2 in. 

kw = thermal conductivity of wire ~ 14 Btu/hr ft °F 
r = radius of wire = 0.0015 in. 

R = radius of assumed insulating cylinder 
= 1 in. for center thermocouple 
= Vi in. for midradius thermocouple 

The thermal conductivity k of the test material was 
taken as 10" 1 Btu/hr ft °F in air and (0.9 to 3) X 10~ 3 
Btu/hr ft °F in vacuum (Tables 2-5); T L was taken as 
the initial temperature of the test, TV 

Solving for the relative temperature change at the mid- 
length produced by the couple, one finds 


atj = T 2 ~ T » = r 2 ~ To 

T 2 - Ti (T 2 - To) - (T 2 - To) 


1 


1 + 


T h - To 
To ~ T 2 


1 


14-2 sinh 2 



(A-10) 


In air, this gives At/ = 2 X 10~ 27 at center and 5 X 10" 29 
at midradius. In vacuum, the corresponding center value 
is 5 X 10~ 3 or 4 X 10~ 5 , depending on whether test ma- 
terial conductivity is taken as 0.9 X 10“ 3 or 3 X 10 -3 
Btu/hr ft °F; the midradius value is 4 X 10~ 3 or 2 X lO* 5 , 
Multiplying by the ratios of — 100 ( Aa/a) /A 17 from Sec- 
tion II-B, percentage errors in diffusivity and conductivity 
are obtained. The errors given in Table A-l correspond 
to the largest values so obtained, and are thus conserva- 
tive. Moreover, the calculation of temperature differences 


due to the thermocouples is based on steady-state con- 
ditions; these are not attained during a test. The transient 
temperature differences will be smaller, so the errors 
quoted are again conservative. 

It may be noted that because sinh x (Eq. A-10) may 
vary rapidly with x, the errors arising from heat flow in 
the thermocouples are very sensitive to experimental 
conditions. For example, if the thermocouple wires had 
been 0.012 instead of 0.003 in. in diameter, the diffusivity 
and conductivity determinations could be in error by a 
factor of 2. 


F. Errors Arising From Heat Flow Through the 
Bottom Closure 

Heat flow through the bottom closure affects the cal- 
culated values of the thermal properties only insofar as it 
introduces axial heat flow through the ends of thermo- 
couples and dust. Errors arising from such axial flow have 
been estimated above. 


G. Errors Arising From Lack of Temperature 
Equalization at Start of Test 

In a series of tests on a given specimen, complete 
equalization of temperature was not always attained 
after one step in copper cylinder temperature and before 
a second step. This would not affect conductivity values 
calculated by digital computer, as the measured radial 
distribution of temperature at the start of the test was 
inserted into the computer. For manual calculations, 
however, uniform temperature at the start of each test 
was assumed. In the worst cases, the initial difference 
between center thermocouple and copper cylinder ther- 
mocouple amounted to about 10 °F, between midradius 
thermocouple and copper cylinder thermocouple to about 
6°F. For 100°F temperature steps, the resulting errors in 
U are about 10% and 6%. Corresponding errors in dif- 
fusivity, when manual computation was used, are given 
in Table A-l. 


H. Errors Arising From Variation of Diffusivity 
with Temperature 

In the manual calculations, it is assumed that the dif- 
fusivity is constant during a test. Results of several tests 
on a sample (Fig. 7-9) show, however, that the diffusivity 
varies perhaps 6% per 100°F in air and 10% in vacuum. 
The calculation gives some sort of average value for the 
temperature range of the test; this is estimated to differ 
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from the diffusivity at the mean temperature by 10 to 
25% of the above variation. 

This error does not arise with machine computation, in 
which conductivity is considered to be a function of 
temperature. 

/. Errors Arising From Variation of Specific Heat 
with Temperature 

An average value of the specific heat over the test 
range was used in manual calculations of conductivity. 
Experimentally, the specific heat varied about 10% per 
100 °F (Fig. 5). The conductivity calculated may differ 
from that at the mean temperature by 10 to 25% of this 


variation. In machine calculations, the specific heat was 
taken as a function of temperature. 

J. Errors Arising From Uncertainties in 
Specific Heat 

The experimental error in measuring specific heat, in 
the range 0 to 200°F, is estimated at 5%. For some of 
the conductivity determinations, specific heat down to 
— 100°F was needed. The extrapolation (Fig. 5) is esti- 
mated to be good to 10% . 

Resulting errors in calculated conductivity are equal 
to those in specific heat. No error in diffusivity is intro- 
duced, as specific heat values were not used in the 
diffusivity determination. 
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